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In this article we examine the hypothesis that the highest energy cosmic rays are complex nuclei 
from extragalactic sources. Under reasonable physical assumptions, we show that the nearby metally 
rich starburst galaxies (M82 and NGC 253) can produce all the events observed above the ankle. 
This requires diffusion of particles below 10'^° eV in extragalactic magnetic fields _B « 15 nG. 
Above 10^® eV, the model predicts the presence of significant fiuxes of medium mass and heavy 
nuclei with small rate of change of composition. Notwithstanding, the most salient feature of the 
starburst-hypothesis is a slight anisotropy induced by iron debris just before the spectrum-cutoff. 
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I. INTRODUCTION 

The cosmic ray (CR) spectrum above 10^" eV (where 
the Sun's magnetic field is no longer a concern) can be 
described by a series of power laws with the flux falling 
about 3 orders of magnitude for each decade increase in 
energy [Q . Above 10^^ eV, the flux becomes so low that 
direct measurements using sophisticated equipment on 
satellites or high altitude balloons are limited in detector 
area and in exposure time. Ground-based experiments 
with large apertures make such a low flux observable after 
a magnification effect in the upper atmosphere: the inci- 
dent cosmic radiation interacts with atomic nuclei of the 
air molecules and produces extensive air showers which 
spread out over large areas. Continuously running moni- 
toring through ingenious installations has raised the max- 
imum observed primary particle's energy to higher than 
1020 eV §. 

While theoretical subtleties surrounding CR acceler- 
ation provide ample material for discussion, the debate 



about the origin of CRs up to the knee (~ 10 



15.5 



3V) 



has reached a consensus that they are produced in su- 
pernova explosions js) . The change of the spectral index 
(from —2.7 to —3.0) near the knee, presumably reflects a 
change in origin and the takeover of another, yet unclear 
type of source. The spectrum steepens further to —3.3 
above r-^ 10^^'^ eV (the dip) and then flattens to an index 
of —2.7 at ~ 10^*'^ eV (the ankle). A very widely held in- 
terpretation of the modulation features is that above the 
ankle a new population of CRs with extragalactic ori- 
gin begins to dominate the more steeply falling Galactic 
population [Q . The origin of the extragalactic channel is 
somewhat mysterious. 

CRs do not travel unhindered through intergalactic 
space, as there are several processes that can degrade the 
particles' energy. In particular, the thermal photon back- 
ground becomes highly blue shifted for ultrarclativistic 
protons. The reaction sequence p^ — > A+ -^ n'^p effec- 
tively degrades the primary proton energy providing a 
strong constraint on the proximity of CR-sources, a phe- 
nomenon known as the Greisen-Zatsepin-Kuz'min (GZK) 
cutoff |l. Speciflcally, fewer than 20% of lO^^-^ eV 
(10^*^ eV) protons can survive a trip of 18 Mpc (60 Mpc) 



M. A heavy nucleus undergoes photodisintegration in 
the microwave and infra-red backgrounds; as a result, 
iron nuclei do not survive fragmentation over compara- 
ble distances M. Ultra high energy gamma rays would 
travel even shorter paths due to pair production on radio 
photons |8|. 

In order to analyze the effect of energy losses in the 
observed spectrum, it is convenient to introduce the ac- 
cumulation factor /ace, deflned as the ratio of energy- 
weighted fluxes for low (lO^^-'^ eV - lO^^-^ eV) and high 
(> 10^" eV) energy CRs above the ankle. In the case 
where the cosmic rays are protons from a uniform dis- 
tribution of sources active over cosmological times, the 
cutoff due to the photopion processes relates the accumu- 
lation factor to a ratio of the GZK distances M and leads 
to /ace ~ 100. A similar value for face is obtained for nu- 
clei due to photodisintegration. Therefore, in the case of 
ordinary baryonic CRs, if Earth is located in a typical 
environment and all CR-sources have smooth emission 
spectra, the observed spectrum above the ankle should 
have an offset in normalization between low and high 
energy given by /ace- To reproduce the recorded spec- 
trum (i.e., /ace ^ 1), the power of nearby sources (say, 
10 Mpc or so) should be comparable to that of all other 
sources (redshift z > 0.5) added together. This condi- 
tion imposes stern constraints on models describing the 
origin of baryonic ultra high energy CRs. For instance, 
"top down" models (with hard injection spectra ex E~^) 
|lO| would fail to reproduce the detected population of 
CRs below the GZK energy by more than an order of 
magnitude. For models that rely on GZK-evading mes- 
sengers mW , /ace depends on the details of the model. In 
the case of messengers which can induce showers across 
the entire energy spectrum, one expects enhancement on 
the low energy side only from the baryonic component, 
and /ace depends on the interaction length in the cos- 
mic microwave background and on the relative energy 
spectra at the source. For messengers whose attenuation 
length is comparable to the horizon, and which do not 
shower at the lower energies, /ace "^ 1- Summing up, 
the smoothness of the observed CR spectrum suggests, 
as the simplest explanation, that nearby sources should 
be signiflcantly more concentrated or more powerful than 



average. 

On a different track, any candidate model address- 
ing the origin of ultra high energy CRs should properly 
match the main features of the observed extensive air 
showers. As the cascade develops in the atmosphere, the 
number of particles in the shower increases until the sec- 
ondary particles' energy is degraded to the point where 
ionization losses dominate, and the density of particles 
starts to decline. The number of particles as a function 
of the amount of atmosphere penetrated by the cascade 
in g cm~^ becomes a smooth curve, the so-called "lon- 
gitudinal profile". The atmospheric depth at which the 
shower reaches its maximum size is referred to as the 
depth of shower maximum Xmax, and is often regarded 
as the most basic parameter of the shower. It increases 
with primary energy as more cascade generations are re- 
quired in the cooling of secondary products. For a given 
total energy X^ax is related to the energy per nuclcon 
of the shower progenitor. Unfortunately, extracting in- 
formation on the nature of the primaries from the air 
shower they produce has proved to be exceedingly diffi- 
cult. The most fundamental drawback is that the first 
few cascade steps are subject to large inherent fluctu- 
ations and consequently this limits the event-by-event 
mass resolution of the experiments. In addition, the cen- 
ter of mass energy of the first interactions is well beyond 
those reached in collider experiments. Thus, one needs 
to rely on hadronic interaction models that attempt to 
extrapolate (using different mixtures of theory and phe- 
nomenology) our understanding of particle physics. 

An analysis of the histogram of Xmax, observed by the 
Fly's Eye experiment, indicates that there is a significant 
fraction of nuclei with charge greater than unity in the 
energy range (lO^^'^ eV - 10"'^^ eV, and somewhat above) 
p2[ . To examine the situation above the GZK-energy, 
simulations of giant air shower evolution have been per- 
formed by means of the code AIRESQ (version 2.1.1) [[l3[ . 
Several sets of protons and iron nuclei were injected at 
100 km above sea level. The geomagnetic field was set 
to reproduce that in the Utah desert. All shower par- 
ticles with energies above the following thresholds were 
tracked: 750 keV for gammas, 900 keV for electrons and 
positrons, 10 MeV for muons, 60 MeV for mesons and 
120 MeV for nucleons and nuclei. The charged multi- 
plicity, essentially electrons and positrons, was used to 
determine the number of charged particles and the loca- 
tion of the shower maximum by means of 4-parameter 
fits to the Gaisser-Hillas function. In Fig. 1 we show 
the evolution of X^ax above the GZK-energy for pro- 
tons and iron nuclei. For comparison, we also show the 
depth of shower maxima of the two highest energy events 
recorded by Fly's Eye and HiRes experiments |lj]. The 
observed values of X^ax are consistent with both proton 
and iron primaries, perhaps (speculatively) suggestive of 
a medium mass nucleus |T5[ . 

In the search for the trans-GZK-sources another ob- 
servable that one has to take into account is the CR ar- 
rival directions. The observed events above « 10^° eV 
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FIG. 1. Simulated depth of shower maxima of two different 
particle species superimpossed over experimental data. 



are distributed widely over the sky, with no plausible 
counterparts (such as sources in the Galactic Plane or 
in the Local Supercluster). Moreover, the data are con- 
sistent with an isotropic distribution of sources in sharp 
contrast to the anisotropic distribution of light within 50 
Mpc. At first glance, this seems to contradict any expla- 
nation based on nearby sources. However, if the highest 
energy CRs are heavy nuclei, one cannot yet rule out 
that extragalactic/galactic magnetic fields could tangle 
up the particle paths, camouflaging the exact location 
of the sources. Intergalactic field strengths and coher- 
ent lengths are not well established, but it is plausible to 
assume that fields have coherent directions on scales of 
^ ~ 0.5 — 1 Mpc. A recent estimate on the extragalac- 
tic magnetic field in the neighborhood of the Milky Way 
suggests i? > 10 nG [M. However, magnetic fields of 
a few /LtG structured in cells of ^ 1 Mpc cannot be ex- 
cluded [|l^. Such a large field would completely defiect 
trans-GZK proton orbits. For a CR nucleus of charge Ze 
in a magnetic field B^q = B/10^^ G, the Larmor radius 
is 
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where Eis = E/10^^ eV. Therefore, the assumption that 
the giant air showers E > 10^° eV were triggered by 
heavy nuclei implies ordered (£ ~ 1 Mpc) extragalactic 
magnetic fields B^g < 15 (at least in the outskirts of 
the Galaxy), or else nuclei would be trapped in magnetic 
subdomains suffering catastrophic spallations [ [l8[ . 

In a previous paper p9| |, there was explored the hy- 
pothesis that CRs above the ankle are (mostly) protons 
from the nearby radio galaxy Centaurus A [^, acceler- 
ated in a "hot spot" in the northern middle lobe. These 
protons were assumed to diffuse in a magnetic field of 
0{idG) over the transit distance ~ 3.4 Mpc to Earth. 
In this work, in view of the possibility raised by the air 
shower profiles, we examine an alternate hypothesis: that 
the composition of CRs above the ankle are largely heavy 
nuclei which originate in two nearby sources, and tra- 
verse extragalactic magnetic fields which conform to the 
restriction outlined in the previous paragraph. As will 
be seen, this hypothesis (as well as the previous one in- 
volving Gen A), will be subject to specific testing in the 
coming array of high statistics cosmic ray observations. 



II. STARBURST-HYPOTHESIS 

If the trans-GZK particles are heavy nuclei, then the 
nearby (~ 3 Mpc |Q) starburst galaxies M82 {I = 
141°, 6 = 41°) and NGC 253 {I = 89°, 6 = -88°) 
would probably be the sources of most ultra high energy 
CRs observed on Earth. Starbursts are galaxies under- 
going a massive and large-scale star formation episode. 
Their characteristic signatures are strong infrared emis- 
sion (originated in the high levels of interstellar extinc- 



tion), a very strong HII- region- type emission-line spec- 
trum (due to a large number of O and B-type stars), 
and a considerable radio emission produced by recent su- 
pernova remnants (SNRs). Typically, the starburst re- 
gion is confined to the central few hundreds of parsecs of 
the galaxy, a region that can be easily 10 or more times 
brighter than the center of normal spiral galaxies. In 
the light of such a concentrated activity, the existence of 
galactic superwinds is not surprising |2^ . 

Galactic-scale superwinds are driven by the collective 
effect of supernovae and massive star winds. The high 
supernovae rate creates a cavity of hot gas (~ 10^ K) 
whose cooling time is much greater than the expansion 
time scale. Since the wind is sufficiently powerful, it can 
blow out the interstellar medium of the galaxy avoid- 
ing it remaining trapped as a hot bubble. As the cavity 
expands a strong shock front is formed on the contact 
surface with the cool interstellar medium. The shock 
velocity can reach several thousands of kilometers per 
second and ions like iron nuclei can be then efficiently 
accelerated in this scenario up to ultra high energies by 
Fermi's mechanism pl| . 

In a first stage, ions are diffusively accelerated at sin- 
gle supernova shock waves within the nuclear region of 
the galaxy. Energies up to ^ 10^^ eV can be achieved 
in this step |g[. Heavy nuclei arc not photodissociated 
in the process despite the large photon energy densities 
(mostly in the far infrared) measured in the central re- 
gion of the starburst. The escape of the CR outfiow is 
convection dominated. In fact, the presence of several 
tens of young SNRs with very high expansion velocities 
and thousands of massive O stars (with stellar winds of 
terminal velocities up to 3000 km s~^) must generate 
collective plasma motions of several thousands of km per 
second. Then, due to the coupling of the magnetic field 
to the hot plasma, the magnetic field is also lifted out- 
wards and forces the CR gas to stream along from the 
starburst region. Most of the nuclei escape in opposite 
directions along the symmetry axis of the system, as the 
total path traveled is substantially shorter than the mean 
free path ||ltl. 

Once the nuclei escape from the central region of the 
galaxy (with energies of ~ 10^^ eV) they are injected 
into the galactic-scale wind and experience further ac- 
celeration at its terminal shock. For this second step in 
the acceleration process, the photon field energy density 
drops to values of the order of the microwave background 
radiation (we are now far from the starburst region) , and 
consequently, iron nuclei are safe from photodissociation 
while energy increases from ~ 10^^ to 10^° eV. In terms 
of parameters that can be determined from observations, 
the nucleus maximum energy is given by |2l|]. 
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where Esw ~ 2.7 x 10*^ erg s ^ is the superwind kinetic 
energy flux and M — 1.2M0 yr^^ is the mass flux gener- 



ated by the starburst ||22l • The age To„ can be estimated 
from numerical models that use theoretical evolutionary 
tracks for individual stars and make sums over the entire 
stellar population at each time in order to produce the 
galaxy luminosity as a function of time |23] . Fitting the 
observational data, these models provide a range of suit- 
able ages for the starburst phase that goes from 50 Myr 
to 160 Myr |g2[. These models must assume a given ini- 
tial mass function (IMF), which usually is taken to be a 
power-law with a variety of slopes. Recent studies have 
shown that the same IMF can account for the properties 
of both NGC 253 and M82 ||. Besides, a region (re- 
ferred to as M82 "B") near the galactic center of M82, 
has been under suspicion to be a fossil starburst site in 
which an intense episode of star formation occurred over 
100 Myr ago |24|. The derived age distribution suggests 
steady, continuing cluster formation at a modest rate at 
early times (> 2 Gyr ago), followed by a concentrated for- 
mation episode 600 Myr ago and more recent suppression 
of cluster formation. In order to get some estimates on 
the maximum energy, let us assume B '^ 50/zG, a choice 
consistent with observation [p5|. Replacing in Eq. (0) 
in favor of these figures, already for Ton = 50 Myr one 
obtains 
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(3) 



Now, we could use the rates at which starbursts in- 
ject mass, metals and energy into superwinds to get an 
estimate on the CR-injection spectra. Let us introduce 
e, the efficiency of ultra high energy CR production by 
the superwind kinetic energy flux. Using equal power per 
decade over the interval 10^*-^ eV < E < 10^°-^ eV, we 
obtain a source CR-luminosity 
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where the subscript "0" refers to quantities at the source. 
The density of CRs at the present time t of energy 
< 10^° eV at a distance r from a source (assumed 
to be continuously emitting at a constant spectral rate 
dNo/dEdt from time ion until the present) is mm 
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where D{E) stands for the diffusion coefficient, x 
ADTon/r'^ = Ton/TD, Ton^t- ton, and 



V^ Ji/x V" 



(6) 



In each "scatter", the diffusion coefficient describes an 
independent angular deviation of particle trajectories 
whose magnitude depends on the Larmor radius. CRs 



with energies E < 10-^^^Mpc -Z^-BnG eV remain trapped 
inside cells of size -^Mpc = -^/(l Mpc), attaining efficient 
diffusion when the wave number of the associated Alfven 
wave is equal to the gyroradius of the particle [g^l- It 
may be plausible to assume a Kolmogorov form for the 
turbulent magnetic field power spectrum, this gives for a 
diffusion coefficient 1271 
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For T„. 



oo, the density approaches its time- 



independent equilibrium value rieq, while for Ton = td 



r^/AD, n/ric 



0.16. To further constrain the param- 



eters of the model, we evaluate the energy-weighted ap- 
proximately isotropic proton flux at 10^^ eV, which lies 
in the center of the flat "low energy" region of the spec- 
trum: 
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where /* = /m82 + ^ngc 253- In the second line of the 
equation we used B^g = 15, ^Mpc = 0.5, and an average 
Z = 20. We fix 



6/^ = 0.013, 



(9) 



after comparing Eq. (^ 



E^JobsiE) = 1024-5 eV2 



to 
-2 



the 

-1 



observed CR-flux: 
m - s"" sr ^ |jl|. Note that the 
contribution of /m82 and /ngc 253 to /* critically depends 
on the age of the starburst. In Fig. 2 we show the re- 
lation "starburst-age/superwind-efficiency" derived from 
Eq. (I). We have assumed that both M82 and NGC 253 
were active for 115 Myr (e « 10%), beyond this epoch 
CR-emission must be associated to M82 "B" . 

Above > lO^^'^ eV iron nuclei do not propagate diffu- 
sively. Moreover, the CR-energies get attenuated by pho- 
todisintegration off the microwave background radiation 
and the intergalactic infrared background photons. The 
disintegration rate of ^^Fe as a function of the Lorentz 
factor r can be parametrized as follows pS], 



i?(r) = 3.25x io-er-"-643 

X exp(-2.15 X 10^°/r) s"^ 

if r e [1. X 10^3.68x lO^''], and 

R{r) = 1.59 X lo^^^r 
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if F e [3.68 X 10^", 1. X 10"]. At this stage it is worth- 
while to point out that the knowledge of ^^Fe effective 
nucleon loss rate is enough to obtain the corresponding 
value of R for any other nuclei |7| , 
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FIG. 2. Age of the starbursts as a function of the efficiency 
of CR-production e. 



This means that A{t) is an exponential function of time 
with an e-folding time of 56 dA/dt\p^ . Now, since the 
emission of nucleons is isotropic in the rest frame of the 
nucleus, the average fractional energy loss is equal to 
the fractional loss in mass number of the nucleus (i.e., F 
is conserved). The relation that determines the energy 
attenuation length of iron as a function of the time flight 
t is then 



E{t) = 938 A{t) r MeV 



Eq exp 



56 



(12) 



where Eq = 938 Aq T MeV, denotes the nucleus' emis- 
sion energy [ p9| . This relation imposes a strong con- 
straint on the location of nucleus-sources: note that less 
than 1% of iron nuclei (or any surviving fragment of their 
spallations) can survive more than 3 x 10^^ s with an en- 
ergy > lO^O-s eV. 

In the non-diffusive regime, the accumulated deflection 
angle 0{E) from the direction of the source, located at a 
distance d, can be estimated assuming that the particles 
make a random walk in the magnetic field [pO| 
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where E20 = E/10'^° eV. Therefore, if B - 15 nG all di- 
rectionality is lost. The resulting time delay with respect 
to linear propagation is given by 
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and the total travel time is 
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As an example, we apply these consideration to the 
highest energy Fly's Eye event. Including statistical and 
systematic uncertainties, the energy of this event is 3.2 ± 
0.9 X 10^" eV. Eqs. ( p^ and ^ relate the uncertainty 
in energy to the uncertainty in the attenuation time: 
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/,From these considerations, we find that the upper limit 
on the transit time for a nuclear candidate for the highest 
energy Fly's Eye event is ~ 6 x 10^** s. The arrival di- 
rection of the highest energy Fly's Eye event is 37° from 
M82 pg. With d ~ 3 Mpc and 6* = 37°, we find from 
Eq. ( p^ a transit time t ~ 3.4 x lO^"* s, well within the 
stated upper limit P^ . 

For average deflections of 60°, the time of flight is 
~ 3.9 X 10^^ s, and consequently there is a sharp end 
of the CR-spectrum near the maximum observed energy. 
Indeed, this leads to a slight anisotropy just before the 
cutoff, and eventually to a north-south asymmetry on the 
tail [Q. It is rather difficult to assess whether events 
with energies > 10^°'^ eV are plausible. This would re- 
quire an event-by-event analysis because the maximum 
energy strongly depends on Tdeiay We do not attempt 
to make yet another estimate in the present article and 
only note that the energy- weighted flux beyond the GZK- 
energy (due to a single M82 flare) p4| 
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is easily consistent with observation Q. The analytical 
study outlined in this paper should be followed up by 
numerical studies, especially to refine our estimates above 
the GZK-energy. 



III. CONCLUSION 

We have shown that the nucleus-emitting-sources of 
our backyard have enough power to produce all CRs ob- 
served above the ankle. Starburst galaxies can accelerate 
iron nuclei above the GZK-energy if a two-step process is 
involved. The crucial point is that for energies > 10^^ eV, 
acceleration occurs in the terminal shock of the starburst 
superwind, well outside the problematic central region. 

Below 10^^ eV, the distribution of the CR arrival di- 
rections is expected to be completely isotropic because 
of Kolmogorov diffusion in ordered (fivipc = 0.5) extra- 
galactic magnetic fields 10 < B^g < 15- In addition, 
(de)magnification of the fluxes by lensing effects are ex- 
pected due to deflections in the regular Galactic magnetic 
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field 1^. Furthermore, medium mass and heavy nuclei 
with energies < 10^^'^ eV would also have an isotropic 
distribution in the sky. 

On the other hand, ultra high energy {E > 10""' e 
light nuclei {Z < 10) do not propagate diffusively. How- 
ever, as can be seen in Eq.(y), light nuclei are not co- 
piously accelerated to ultra high energies in the star- 
burst. Note also that the nucleons (a-particles) emit- 
ted in the photodisintegration process have energies well 
below 10^^ eV and consequently do not produce any 
anisotropy. 

It should be noted that several groups Q have re- 
cently reported evidence of clustering (6 doublets and 1 
triplet, with the chance probability 7 x lO^"' I23l)- Mag- 
netic focussing in the magnetic field structure could - in 
principle - account for directional clustering to explain 
the current sparse data [Q. However, if not a statisti- 
cal fluctuation, and clusters are well established in very 
much larger data sets, they would constitute a serious 
objection to the model. 

All in all, within this scenario almost all CRs above 
the GZK-energy would be medium mass and heavy nu- 
clei, yielding a small rate of change of composition (in 
agreement with observation |p9|). In addition, the model 
predicts a slight anisotropy above lO^^** eV produced 
by non-diffusive iron debris. This makes the model ca- 
pable of being proven false. The limited statistics in 
the observed data make it impossible to definitively test 
the "starburst-hypothesis" at this time. The coming 
avalanche of high quality CR-observations ||l[ premisses 
to give the final verdict on these speculations. 
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